Electronic structure of edge and vortex states in chiral mesoscopic superconductor. 
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We study a subgap quasiparticle spectrum in a mesoscopic disk of chiral superconductor. We 
find an exact expression for the spectrum of surface states localized at the disk edge. Considering 
an Abrikosov vortex placed at the center of a superconducting disk we investigate the spectrum 
transformation near the intersection points of surface and vortex anomalous energy branches . The 
resulting splitting of the anomalous branches due to the hybridization of edge and vortex states is 
determined by an external magnetic field and can lead in particular to the formation of a set of 
minigaps in the quasiparticle spectrum. Tuning the external magnetic field makes it possible to 
control the width of minigaps and the positions of corresponding density of states singularities at 
the minigaps edges. 

PACS numbers: 74.25.-q, 74.78.Na 



1. Introduction. Recently, a considerable attention 
has been devoted to the investigation of a chiral super- 
conducting state which is proposed to realize in p-wave 
triplet superconductor Sr2Ru0 4 (see Ref. [jij and ref- 
erences therein). The chiral superconductivity can be 
associated with a formation of Cooper pairs with a non- 
zero orbital angular momentum. In this case a value of 
chirality x is determined as a projection of the angular 
momentum of Cooper pairs on the z axis. Generally, the 
superconducting order parameter is triplet (singlet) for 
odd (even) chirality x anc i is given by A = A e IX F a z 
and A = Aoe lxSp correspondingly where Ao is a bulk 
superconducting gap and a z is a conventional spin opera- 
tor. For x 7^ the phase of the order parameter depends 
on the direction of the electron momentum in xy plane: 
p = p(cos9 p , sin# p ). The important consequence of this 
fact is an existence of surface Andreev bound states[l,|3]. 
They appear in the vicinity of scattering interfaces be- 
tween a superconductor and an insulator, if the order 
parameter phase takes different values for the incident 
and reflected quasiparticles (QP) with different momen- 
tum directions. The formation of Andreev bound states 
increases the local density of states at the surface of a su- 
perconductor resulting in zero-bias a conductance peak 
anomaly observed in tunneling spectroscopy of high-T c 
cuprates with a d-wave symmetry of superconducting 
pairing [3] as well as in p-wave triplet superconductor 
Sr 2 Ru0 4 @. 

Under an applied magnetic field, generating screening 
current and vortices the spectrum of surface states ac- 
quires a Doppler shift, leading to the splitting of the zero- 
bias conductance peak Recently in work 8] it was 
proposed that the Doppler shift effect should lead to the 
chirality selective influence of the magnetic field on the 
surface states in chiral p-wave superconductor (\x\ = 1), 
such as Sr2RuOi. The QP density of states (DOS) near 
the flat surface was shown to depend on the orientation 
of the magnetic field with respect to the z axis as well 
as on the vorticity in case when an Abrikosov vortex is 



pinned near the surface of a superconductor. 

If an Abrikosov vortex is situated close to the boundary 
of a superconductor, then in addition to the Doppler shift 
effect 8] it is necessary to take into account the hybridiza- 
tion of edge modes and low-energy vortex core states [3] . 
In this case spectrum modification is determined by the 
overlapping of QP wave functions localized near the sur- 
face and near the vortex core. The characteristic local- 
ization length of subgap QP wave functions is determined 
by a superconducting coherence length £. Therefore the 
hybridization of vortex and edge states should be partic- 
ularly important in mesoscopic superconducting samples 
of the size of several £. 

Let us consider a model problem when a supercon- 
ducting sample has ideal disk geometry in xy plane. In 
this case the spectrum of edge states can be expressed in 
terms of the angular momentum fj, which is conserved due 
to the axial symmetry. We assume the validity of a qua- 
siclassical approach so that one can consider QP motion 
along trajectories, i.e. straight lines along the direction of 
QP momentum p = p(cos9 p ,sm9 p ). Employing analogy 
with a point Josephson junction the spectrum of surface 
states can be written as follows: E s = — Ao cos(Ay>/2), 
where < Aip < 2ir is the difference between gap func- 
tion phases seen by the incident and reflected QPs. Un- 
der the reflection of a trajectory at the disk boundary 



the angle 9 P transforms as 9 p — > 9 p 



2 arcsin(6/_R), 



where b = —fi/kp is a continuous impact parameter, i.e. 
a distance from the trajectory to the disk center, Uf is a 
Fermi wave number and R is a disk radius. In case when 
there is no circulating superconducting currents we ob- 
tain A<p = [x{tt + 2 arcsin(6/i?))] mod (2tt) yielding a 
set of anomalous energy branches corresponding to the 
edge states (13, EH: 



E 



(1) 



where j = l...\x\ and fj,j — (kpR) sm(7rnj/2x). The inte- 
ger index rij from the interval — \x\ < n j < Ixl is chosen 
so that the combination x ~ n j to be odd. As noted in 
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Ref. [lj| the spectrum of edge states (JTJ) is analogous to 
the general spectrum of quasiparticles localized within 
a vortex core [12]. The inter level spacing for a particu- 
lar anomalous branch u) s j = xA /(fcpi? cos(nnj /2x)) is 
much smaller than the bulk superconducting gap Ao pro- 
vided kpR ^S> 1, therefore the anomalous branches can 
be considered as functions of a continuous impact param- 
eter b = —fi/kp- In case of even chirality \ all energy 
branches cross the Fermi level at finite impact parame- 
ters b = —/ij/kF, and for the odd \ there exists an energy 
branch with fij = 0, crossing the Fermi level at b = 0. 

If an Abrikosov vortex is placed at the center of a su- 
perconducting disk there appears another anomalous en- 
ergy branch associated with the spectrum of vortex states 

E v (fi) = -[iw v . (2) 

Here to v ~ r]Ao/(kF£,)i where £ = HVf / Aq is a supercon- 
ducting coherence length and Vf = frkF/m is a Fermi 
velocity. The value of vorticity rj = ±1 is determined by 
a sign of the superfluid velocity circulation in the counter- 
clockwise direction around the vortex core. The angular 
momentum /i in Eq.([2]) is integer (half-integer) for the 
odd (even) chirality value 0]. 

Considered as continuous functions of a quasiclassical 
impact parameter b = —fi/kp the spectrum branches 
E v (b) and E S j(b) intersect at the certain points b = bj. 
The splitting of energy levels at the degeneracy point oc- 
curs due to the hybridization of vortex and edge states 
and can be estimated using the perturbation method 
for an almost degenerate two-level system (see Ref. [Hi]), 
which yields the secular equation: 

[E — E s j(b)][E - E v (b)] = J 2 , (3) 

where the factor J is determined by the overlapping of the 
corresponding wave functions. Using a Taylor expansion 
E v (b) = E v (bj) + E' v (b 3 ){b - bj) and E sj (b) = E sj (bj) + 
E' s j(bj)(b — bj) one can see that the scenario of branch 
splitting depends on the slopes of energy branches E' v — 
dE v /db and E' s j = dE S j/db at the intersection point b = 
bj . In case when the signs of the slopes are opposite there 
appears a minigap in the QP spectrum. The minigap 
width i.e. the minimal energy spacing between QP levels 
corresponding to the different energy branches can be 
found from Eq.Q as follows: 

5E = 2\J\ ^\E> v E> sj \/{\E> v \ + \E' s j\). (4) 

Otherwise, when E' v and E' g j have the same sign, the mini- 
gap width is always equal to zero. Generally the energy 
branches E v {b) and E S j(b) cross at \bj\ ~ £ and \E\ ~ A , 
where the spectra are not described by the expressions 
(|1I2|) . For the sake of simplicity further we focus on a par- 
ticular case of chiral p-wave superconductor with \x\ = 1. 
Then there is only one anomalous surface energy branch 



E s (b), crossing the Fermi level at b — simultaneously 
with the vortex energy branch E v (b) . Therefore the spec- 
trum transformation takes place within a domain of small 
energies \E\ <C A , where the surface and vortex states 
are well localized and the overlapping factor can be eval- 
uated as follows | J\ ~ A e~ R ^. In this case the splitting 
of energy branches is shown schematically in Fig.JT]) for 
the opposite and equal signs of chirality \ an d vorticity 
r). 
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FIG. 1: FIG. 1 Shown by solid lines is the spectrum transfor- 
mation caused by the interaction of vortex and surface spec- 
trum branches: (a) chirality and vorticity have different signs 
(v — 1) X = — l)i (b) chirality and vorticity have the same 
signs (rj = 1, x = !)• The dash and dash-dotted curves 
correspond to the non-interacting vortex and surface energy 
branches. 

Following Ref. Q one could expect the surface en- 
ergy branch E s (b) to be modified due to a Doppler shift 
induced by a superfluid velocity circulating around the 
vortex core. However, it is not the case if \x\ = 1 when 
the Doppler shift is totally compensated by an additional 
vortex-induced difference of the order parameter phase 
for the incident and reflected QP. Indeed, in the presence 
of a vortex-induced superfluid velocity v s the Doppler 
shifted spectrum of surface states is given by: E s — 
— A cos(Aiyj/2) + pf • v s , where Pf • v s = —rjA (b/R) 
is a Doppler shift energy at the disk edge and the phase 
difference is Aip « n + 2(x + rj)b/R. It is easy to see 
that for \x\ = 1 the surface energy branch is still given 
by Eq.|T|) with rij = 0. 

2. Model. In order to investigate in detail the effects 
described above we proceed with a quantitative analysis 
of the QP spectrum in chiral superconductor on the basis 
of Bogoulubov - de Gennes theory: 

where Hq — f 3 [(iW + ef 3 A/c) 2 — p F ]/2m, ^ = (u,v), u 
and v are the amplitudes of the electron and hole compo- 
nents, fi are the Pauli matrices in a particle-hole space, 
A is a gap operator: 

A = {A(r), e ^}, (6) 
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where r is a coordinate operator, A(r) describes the 
spatial dependence of the gap function and {A, B} = 
AB + BA is an anticommutator which provides the gauge 
invariance of A. Here we omit the spin-dependent part of 
the gap operator A, neglecting the spin-orbit interaction. 
Also we neglect the dispersion of QP energy in the direc- 
tion perpendicular to the anisotropy plane xy, assuming 
a cylindrical Fermi surface. The magnetic field is directed 
along the z axis H = —Hzq and for extreme type-II su- 
perconductors we can consider the magnetic field to be 
homogeneous on the spatial scale R and take the gauge 
A = [H, r]/2. Within the quasiclassical approach the 
wave function in the momentum representation can be 
taken in the form: 

+00 

*(P) = 7^ / dse- i(|p| - fifcF)fl/ V(s>0p) , (?) 
k F J 

— 00 

where s is a coordinate along a trajectory. In (s,9 p ) 
representation the expression for the coordinate operator 
in Eq.((6|) has the following form: 



f = sk F /k F + {[k F ,z ],jj} /(2k 2 F ), 



(8) 



where /t = —id/d9 p is an angular momentum operator. 
The equation for ip(s, 9 P ) along a quasiclassical trajectory 
reads: 



d 

tT 3 hV F —ip 

OS 











+ 0) 



where u>jj — \e\H/mc is a cyclotron frequency. The 
terms quadratic in H were neglected in © because we 
consider the distances much smaller than the cyclotron 
radius rjj = V f /ujh- The wave function in the real 
space is expressed from Eq.© in the following way (see 
Refs. 



■in 



de p e lkFrcos(e - e ^ij{rcos(e-9 p ),9 p ). (10) 



For an ideal disk it is convenient to use a polar coordinate 
system (r, 9) with the origin at the disk center. Then, the 
boundary condition at the surface of a superconducting 
disk of the radius R reads: 



V(R,0) = 0. 
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3. Spectrum of edge states. At first we investigate 
the spectrum of edge states solving Eq. © with spatially 
homogeneous order parameter distribution: A(r) = Ao 
and applying the boundary condition (fTTI) . Due to the 
axial symmetry of the superconducting sample we sepa- 
rate the 9 P and s variables in Eq.©: 



(12) 



where fi = n + x/2 is an angular momentum, and n is 
integer. The function G M satisfies the following equation: 

-if 3 hV F -^G fl + A fiG IM = EG tt , (13) 

where E = E + [1 (Hu>h/2). In order to apply the bound- 
ary conditions (TPTj) we evaluate the integral in (JTDJ) using 
the stationary phase method. For a given /1 we obtain: 



_ i(fcj?s*-7r/4) 



V;( S *,0i) + e^/ 4 - fc ^*V(- S *,0 2 ), 



where s* = \J R 2 — (/i/k F ) 2 and the stationary phase 
points are given by 9\ = 9 + arcsin(/i/fci?i?) and 92 = 
9 + 7r — arcsin(/i//cpi?). Thus, we obtain the boundary 
condition for the function G^s): 



G M (s*) = e iQ -^G M (-s*), 



(14) 



where a — — 2 avcsm(p,/k F R)] — 2k F s* — n/2 and 
tp = x[arcsin(/z/fci?i?) — ir/2]. 

The general solution of Eq. (fl~3|) can be written as fol- 
lows: 



1 

e i 7 



1 



(15) 



where c, d are the scalar coefficients, 7 = arccos(i?/Ao) 
and q = ^/Aq — E 2 / Aq. Then, from the boundary con- 
dition (fT4|) we obtain the expression for the spectrum of 
edge states: 



E = A M/V1 + M 2 - fi(huj H /2) 



(16) 



where M — coth(2<7s*/^) cotcp — cosa/ sinh(2gs*/^). It is 
possible to evaluate Eq. (Til))) to find an explicit expression 
for the energy levels lying much lower than the supercon- 
ducting gap. For simplicity we start our analysis with the 
case of a zero magnetic field. Considering the low ener- 
gies \E\ -c Ao we obtain that the spectrum consists of 
|x| energy branches: 



E S j(lj) = -(/x - fij)w s j + A 



(— l) kj cos a 
sinh(2s*/0 ' 



(17) 



where cu S j — x^o coth(2s* /^)/(k F Rcos{-Krij /2\)) and 
fj,j = (k F R) sin(7rnj/2x). The integer index rij from the 
interval — |x| < fij < \x\ is chosen so that the combi- 



nation x 



to be odd: x 



2/c-j 



1 . From 



Eq. (fT7j) one can see that the energy levels are the oscillat- 
ing functions of a disk radius with a period 8R = ir/k F . 
The amplitude of energy levels oscillations is larger than 
the interlevel spacing provided the disk radius is smaller 
than the critical value R c determined by the condition 
uj S j ~ Ao/ sinh(2s*/£). For the typical values of the pa- 
rameter ~ 10 2 — 10 3 we obtain R c /£, ~ 3 — 5. Note 
that in case |x| = 1 Eq. pT|l is completely analogous to 
the expression obtained in Ref. [1 61] for the spectrum of 
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vortex core states modified by the normal reflection of 
QP at the surface of s-wave mesoscopic superconductor. 
At R ^> R c the exponentially small oscillating term in 
Eq.(fT7]) can be omitted and we obtain the low-energy 
spectrum of surface states consistsing of a set of anoma- 
lous branches JT]), similar to the spectrum of multiquan- 
tum vortex with a vorticity equal to \ \ 12J . Then, Eq. (JTTj) 
describes an appearance of an energy band due to the in- 
teraction of surface states localized at the opposite ends 
of a trajectory s = ±s*. The bandwidth Aq/ sinh(2s*/£) 
is proportional to the overlap of decaying wave functions 
of edge states. 

Applying a magnetic field H along the z axis one in- 
troduces the shift of surface energy levels /x(fio;///2) in 
Eq. (fPTj) . In fact it is a Doppler shift effect due to the 
Meissner current flowing along the circumference of a 
superconducting disk. The same effect was studied in 
Ref. [§| for the flat geometry of the superconducting sam- 
ple boundary. Let us consider the expression (|17|) in 
more detail for the case of |x| = 1. Taking into ac- 
count the finite external magnetic field we obtain that 
the spectrum is given by Eq. (|17p with uj s \ — lj s = 
xAq / (kpR) + hcuH /2 and the phase of energy oscillations: 
a(ji) = fj,ir-2k F R-ir/2. Note that a(n+2) = a(^) + 27r, 
therefore the spacing between levels corresponding to the 
angular momentum values [i and /i + 2 is 2u> 8 and can be 
neglected within the quasiclassical consideration. On the 
other hand, the spacing between levels corresponding to 
[i and [I + 1 is Aoe~ 2R ^ sxn(2k F R), which can be much 
larger than uj s if R < R c . Thus one can consider two 
continuous branches corresponding to the odd and even 
values of n = /i — 1/2: 

E s {b) = u s k F b ± 2 X A e- 2R/i sin(2fc F i?), (18) 

where b — —fj,/k F is an impact parameter. At R > £, 
the Doppler shift energy k F b{fou)H /2) in Eq. (fl"8]) can sub- 
stantially change the slope of the branches E s (b). In- 
deed, huj H ~ (H/H c2 )A /(k F £), where H c2 ~ </>o/£ 2 is 
the upper critical field and <po = irhc/e is the magnetic 
flux quantum. Therefore the magnetic field of the mag- 
nitude | H | > (£/R)H C 2 can reverse the slope of E s (b). 
Particularly at H — —2x4>o/{£,R) we obtain a dispersion- 
less energy branches E s — ±2Aoe~ 2R ^ cos(2k F R). In 
the Fig.([2|) we show one of the spectrum branches (for 
/i = 2n + 1/2) given by Eq. (fTB|) for the different values 
of magnetic field. Considering the measurable character- 
istics of the QP spectrum we obtain that the density of 
states (DOS) at the Fermi level v(Q) = 1/\lo s \ can be con- 
trolled tuning the magnitude and direction of an external 
magnetic field. 

Analyzing the influence of external magnetic field on 
the surface states spectrum we have neglected a magnetic 
field H s generated by the current, carried by the surface 
states. The density of this current is of the order of the 
critical one for depairing , [ll[ and it flows within the 
surface shell of the width £. Evaluating the magnetic 




FIG. 2: FIG. 1 Spectrum of surface states for \ — — 1 and 
fi — 2n + 1/2. Curves from bottom to top correspond to the 
different values of magnetic field from H = to H = 2<j>o / (£R) 
(corresponding to ui a =0). We choose R = 4£ and &f£ = 200. 

field we obtain H s ~ (£/A) 2 -ff C 2, where A is a London 
penetration length. Thus for extreme type-II supercon- 
ductors (£/A <C 1) the field of the surface current can be 
neglected. 

The considered model with a spatially homogeneous 
gap function A(r) = Ao is adequate only for not very 
large applied magnetic field. Generally, it does not work 
when the field is large enough to suppress a surface bar- 
rier preventing vortex entry H c ~ H C 2(£,/R) [l7|. Cer- 
tainly, in our case a criterion for the vortex formation 
should be sensitive to the orientation of magnetic field 
with respect to the z axis. Now we proceed with the 
analysis of the QP spectrum assuming that a vortex has 
already entered the sample and is placed at the center of 
a superconducting disk. 

4. Vortex induced spectrum transformation. In 
the vicinity of a vortex core the gap function has the 
following form: A(r) = AoD v (r)e l7ie , where rj — ±1 is 
the vorticity and D v (r) is a dimensionless vortex core 
profile. The gap operator A in (s, 6 p ) representation has 
the following form: 

A = A ^jM {(a + w/k F ), e *Cx-H^} . (19) 

To be specific, we choose a model vortex core profile: 
D v {r) = r/^Jr 2 + £ 2 . Once again we can separate the s 
and Op variables: 

ip(s, Op) = e^+'^+^^G^s), (20) 

where /i is integer. We will consider the trajectories 
passing close to the vortex core with impact parameters 
b £, corresponding to \fi\ <§; kp£. Then, from the 
boundary condition (TTTj) we obtain Eq. (fT4")) with s* = R, 
a = (Jtw- 2k F R-n/2 and <p = (\ + i])[n/(k F R) - vr/2]. 
The function G^(s) we satisfies the following equation: 
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d 



'd 



ihV F f 3 i-G „+U(s)G„ = [E — /-W(s) ) G„. (21) 



The matrices U,W(s) are defined at —R < s < R as 
follows: 



U(s) = A D v {s) — fi 



W(s) = A A, (5)7^2 

We assume that the size of the disk is rather large: R > £ 
therefore Eq. (|2"Tj) together with the boundary condition 
(TTJJ describe the interaction of vortex states localized 
near the vortex center i.e. s = and the edge states 
localized at s — ±R. Thus, it is natural to use the tight- 
binding approximation of the wave function to calculate 
the spectrum. To this end we consider a continuation of 
the function G M (s) and coefficients U,W(s) of Eq. (j2Tj) 
to the whole axis — 00 < s < 00 using the periodicity 
conditions: 

G fl (s + 2R) = e ta - lf3ip G^s) 

U, W(s + 2R) = e~ 2lf ^V, W(s). 

We find the solution of Eq. (j2Tj) as a superposition of 
the functions localized at s n = 2nR corresponding to the 
vortex states and localized at d n = (2n+l)R correspond- 
ing to the surface states: 



(22) 



where C\ , C 2 are the arbitrary coefficients and the generic 
terms are: 

V n (s) = e in{a - f3v) e~ K ^ s ~ s ^e"^ /A , 



K v (s-s n )= / D v (s - s n )fl(s - s n ) 



ds 
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and 



S n (s) = e in ( a - f 3f) e -K 3 (s-d n ) e -if 3 7r/4 



1 ' 



f s ds 
K s (s-s n )= D v (s - d n )d(s - d„) — 

where we have introduced the step function $(s) = s/\s\. 

Following the standard tight binding method we sub- 
stitute the solution in the form (|2"2"|) into the Eq. (|2"Tj) . 
multiply by V* and S* from the left and integrate over s, 



taking into account the overlapping of the nearest neigh- 
bor functions. We omit here the details of calculation of 
the corresponding integrals which yields a linear system 
of equations for the coefficients C\ , 

[E-E v {n)]C 1 = -{i/2)e- R 't{l- 



*)C 2 



[E - E.{ji)\ C 2 = ie- H ^(l - e ia )C x . (23) 

The solvability condition for this system is given by the 
Eq.© with J = ^/2A Q er R /^ sin(a/2). The vortex energy 
branch E v (fi) is given by (|2|) with ui v » 0.84^Ao/(fcFC) + 
Hujh/2 and the surface energy branch E s (/i) is given by 
((T|) with rij — and oj s w x^-ol '{kpR) + fkun/2. 

According to the arguments presented in introduction 
in case when the vorticity r\ and chirality x are of the 
opposite signs [Fig.flja], there exists a minigap at the 
Fermi level. In this case the DOS has van Hove singular- 
ities at the minigap edges, i.e. at E = ±5E. Applying 
the magnetic field H along the z axis one can change the 
slope of the intersecting branches E v (b) and E s (b) and 
therefore shift the positions of van Hove singularities ac- 
cording to the Eq.(j4]). For large enough magnetic fields 
\H\ > 2<fro/(£R) the slope of energy branch E s (b) can be 
reversed. In this case as well as in case of the equal signs 
of vorticity r\ and chirality \ the spectrum is gapless as 
shown on the [Fig.{T])b]. 

5. Effect of surface roughness. In conclusion we 
note that in previous sections we analyzed the electronic 
spectrum for a sample with the perfect surface. One 
can expect that the surface roughness would break the 
interference of QP waves and as a result the DOS sin- 
gularity at the minigap edge would be smeared. How- 
ever, it is not the case for a wide class of surface im- 
perfections. Particularly, let us assume that the surface 
of the sample is described by the equation r = R(9), 
where R{6) is a smooth function fluctuating over an 
average value r = Rq, so that \R(6) — Rq\ <C £ and 
\dR/d9\ <g; £. Then, following a procedure developed 
in Ref.[l4j one can find the solution of Eq.((9]) in the 
form of expansion (f2"2"|) with the angle-dependent coef- 
ficients Ci,C 2 (9). Using a tight binding method we ar- 
rive at the system of differential equations for the func- 
tions C\ , C 2 (6) , which coincides with Eq. (f23|) if one takes 
fi = -id/36 and e ia = -i e -^F(R{e)+R(e+Tr)) e i^[i _ Con _ 

sidering the low-energy limit \E\ <C A , it is natural to 
assume that the functions C± <2 (6) consist of a limiting 
number of the lowest angular harmonics with \fi\ <C 
Thus, the coefficients of Eq. ([2"3|) can be averaged over the 
small angular interval to exclude higher angular harmon- 
ics of the rapidly oscillating exponent e - lfc HA(0)+i?.(0+7r)). 
It is important that the right hand sides of equations 
in system (|23[) can not vanish after this averaging, in 
contrast to the analogous problem for the case of s- 
wave superconductor [14j. The resulting spectrum should 
have the form ([3|) with J = (3Aoe~ Ro ^ , where a factor 
P r-j 1 depends on the particular realization of the surface 
roughness. 
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